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SUMMARY 


In this paper the maxima of amplitudes of thickness and loading noise har- 
monics are established when the radial distribution of blade chord, thickness 
ratio, and lift coefficient is specified. It is first shown that only airfoils 
with thickness distribution and chordwlse loading distributions which are 
symmetric with respect to midchord need be considered for finding the absolute 
maxima of thickness and loading noise. The resulting chordwlse thickness and 
load distributions for these maximum noise conditions require infinite slope at 
some points along the chord but otherwise are uniform. It is shown that sweep- 
ing the blades reduces the thickness and loading noise, but there is no optimum 
sweep which generates the lowest noise. 


INTRODUCTION 


In the design of a high tip-speed rotating blade such as a helicopter rotor 
or a propeller, one important acoustic question is: given radial (spanwise) 

load distribution, thickness ratio, and chord distribution of the blade, can 
the maximum of the level of each of the sound harmonics be established? These 
maxima, of course, correspond to the worst possible acoustic design. If these 
maximum levels are kept within acceptable limits, then neither the chordwlse 
load distribution nor the airfoil shape would be of concern in the acoustic 
design. Another question which comes to mind next is whether sweeping the 
blade tips appropriately can result in the lowest possible noise. In this 
paper both of the above questions are studied and answered. 

The starting point of our analysis is the following equation. Let S be 
the mean surface of the blade, h and p be the local thickness and load distri- 
bution on this surface, respectively. Then, the acoustic pressure p' (x,t) in 
the far field is given by the equation 


+Work supported by U.S. Army Research Office (Durham). 
*Work supported under NASA Research Grant NSG 1474. 


373 


(1-a) 


47rr^p'(x,t) = + 


'i'(x,t) = / [pp^] dS 

S ' r ret 


L^(x,t) = 


Pn 

- / ds 

S r ret 


(1-b) 


(1-c) 


The first term in eq. (1-a) is called the thickness noise. This formulation of 
thickness noise was derived by Hanson (ref. 1) and by Farassat (ref. 2), using a 
different approach. The second term in eq. (1-a) is the loading noise. 

It is assumed that the blade system, lying in y^y 2 -plane, is not in motion 

as a whole. That is, only hovering rotors and static propellers are considered 
here. It is also assumed that unsteady loading noise is negligible. For high- 
speed rotating blades, this assumption is justified for observer positions 
where the sources on the blades appear noncompact. Under this assumption, the 
acoustic pressure will be periodic with fundamental frequency based on blade 
passage frequency. For simplicity one blade is considered in the analysis. 

The nth Fourier component of the noise, p'(x), is found from eq. (1-a) to 
be 


47rr p ' (x) 
o*^n 


2 2 y, p,. 

p n 0 ) Y (x) 
o n 


A 

- ikn r. L. (x) 
1 in 


( 2 ) 


In the following analysis, the surface integrals with respect to S, used in 

evaluation of 'i' (- 5 ) and L. (?) , are written in an unconventional manner in chord- 
n in ’ 

wise direction. Written in this form, the effect of sweep can be introduced 
easily. The bounds are obtained in two stages as follows. First, it is shown 
that if the airfoil shape is deformed in such a way that the chordwise distance 
between the points of equal thickness on the airfoil surface is not changed, 
then the thickness noise is maximum if the airfoil is made symmetric with 
respect to a radius of the blade disc at each radial position . A similar result 
holds for the chordwise loading distribution. Therefore to obtain the absolute 
maxima of thickness and loading noise, only airfoil shapes and chordwise load- 
ing with midchord symmetry should be considered. These and the related result 
concerning sweep do not apply at high frequencies due to the mathematical 
limitations of some of the inequalities used in their derivation. The range of 
applicability of these results is, however, wide, particularly in the case of 
helicopter rotors and conventional propellers. 
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In all the examples 
chord of 0.4 m Is used. 
0.95. 


In this paper a rotor blade 
The number of blades Is two 


of 5-m radius and uniform 
and the tip Mach number Is 
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b 

b' 
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g(y).g(y') 
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I, I 
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J 
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SYMBOLS 

function of n (used In eq. (12)) 

function of n (used In eq. (14)), number of blades 
chordwlse variable (see fig. 1), m 
mean chord In blade tip region, m 
blade chord (function of ri) , m 

value of b£b^ where sln(nb/2n) achieves Its peak, m 

coefficient used In eq. (5) 
speed of sound, m/sec 

section lift coefficient (a function of r|) 
arbitrary positive functions 

blade section thickness variable (see fig. 1), m 
maximum thickness of blade section (function of n) » 

Integrals defined In eq. (7) 

Bessel function of first kind on nth order 


k 


wave number, u/c 


1’ In 



n 


P 

^m 

P' 


surface Integral used In calculation of loading noise, and Its 
amplitude of the nth Fourier component of L. (see eq. (1-c) 
and (2)), 1=1, 2, 3 

Mach number along radiation direction 
tip Mach number of the blade 
harmonic number 

unit normal to surface S, direction from pressure to suction 
side of the blade. 1=1, 2, 3 

load distribution of the blade. Pa 
peak section load (function of rj). Pa 
acoustic pressure. Pa 

amplitude of the nth Fourier component of p'. Pa 
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function defined in eq. (10-a) 

|x-y[ , m 

observer distance from center of rotation, m 
(x^-y^)/r, radiation vector 

blade outer and inner radius, respectively, m 

mean surface of the blade 
observer time, sec 

thickness ratio of the blade (a function of n) 
function defined in eq. (10-b) 

period of the sound, sec 

observer position vector, origin at rotation center 

source position vector, origin at rotation center 
geometric angle of attack (function of n) , deg 
azimuthal angle, rad 

(g^+3^)/2 azimuthal angle of point C midway points A and B in 
fig. 1, rad 

functions of q and h or p indicating azimuthal angles of points 
A and B in fig. 1, rad 

arbitrary function (see eq. 8) 

Dirac delta fimction 

angle between axis of rotation and x, rad 

radial position variable, m 

variable defining the degree of blade sweep 

3 

density of the undisturbed medium, kg/m 
source time, sec 

surface integral used in calculation of thickness noise (see 
eq. (1-b)) 

amplitude of the nth Fourier component of 
angular velocity of the blade, rad/sec 
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DERIVATION OF THE BOUNDS 


In this section, our attention will be focused on The manipulations 

for L. are identical to those of 'F . If T is the period of the sound, then 
in n 





V (x, t)e 


inoit 


dt 


(3) 


For our purpose, the volume T(x,t) vri.ll be written as 


R 




'i'(x,t) = ; ndn / dh / [|l-M^l] 


-1 


ds 


R. 


ret 


(4) 


where ri is the radial position and h is the maximum thickness of the airfoil. 

m 

The azimuthal angle is denoted by g. The angles g and g. are the azimuthal 

t ^ 

angles of points A and B, respectively, in fig. (1). Note that 

g^=g^(r|,h), and (i 1 ,h). The only dependence on x and t in eq. (4) comes 

through the integrand so that the time integral in eq. (3) commutes with all 
the integrals in eq. (4). 

We now introduce the source time t in a manner used by Hawkings and Lowson 
(ref. 3). Since dt= [ | 1-M^ | dt , the Doppler singularity in eq. (4) is can- 
celled. Writing t=T+r/c and using the well-known integration with respect to 
T which results in a Bessel function of first kind (ref. 3), we get 


R 


'y^(x) = C / 


R. 

1 


h 

m 

qdn f 

0 


dh f 


-ing 
e J 

n 


(nknslne)dg 


(5) 


Iknr 

where C=(i) e ° is a constant. Let g =(g +g„)/2 and use g„-g =b/n where 

c t Jo Jo t 

b=b(q,h) is the chordwlse distance between points A and B in fig. (1). Then 
Integrating eq. (5) with respect to g results in 


T (x)=-21C f 
" R. 


f nJ (nknslne) 
0 " 


-m3 

c 

e 


• /Tib V 

srn(^) 


dhdp 


( 6 ) 


We note that the angle g^=g^(n,h). We will show that the maximum of 
when only g^ varies corresponds to g^=constant, that is at each radial station 
and for all 0_^h<h^, points A and B should be located S 3 anmetrlcally with respect 
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to one and the same radius of the blade disk. We need the following result. 
If g(^)_^0 and yi'f) is an arbitrary function, then in any region D 

I = |/ e^”'^ g dy| £ / g dy = I (7) 

^ D D 

To prove this, we note that 

- |l 1^ = / / 2 sin^ j [Y(y)-Y(y') ]g(y)g(y')dydy’^ 0 (8) 

' D D 


In fact, strict Inequality holds in most cases since the integrand in eq. 

(8) has to be non-zero only in a small region in DxD. Note that if Y(y)= 
constant, I I^ I = I. 


In eq. (6), J^(nknsine)^0 even for transonic and low supersonic speeds, 

since nknsine<nM^ where is the tip Mach number. If sin(^)>0, then we can 
— t t 2q — 

apply the result proved above. This would require nh/2r\^ or n^2T:ri/b. Since 
for high-speed blades, the tip region is responsible for the generation of the 
noise, a reasonable value for n is n^2Trx- 7R/b'=4.4R/b' where b' is the mean 
chord in the blade tip region. For blades with B blades; we must have 
nB^4.4R/b'. For two-bladed helicopter rotor blades the following result typi- 
cally holds up to twenty fifth harmonics of the blade passage frequency. Apply- 
ing eq. (7) to eq. (6), we get 


I'i (x) < 2 f 
n — 


m 


R. 


qJ (nkqslne) sln(^) dhdq 
n /q 


(9) 


In exactly similar fashion, we can show that 

„ R ^m 

Iq (x)|= It^L- +r»L„ 1< / / J (nkqsine)sin(^)sina(q)dpdq , (10-a) 

n 'ilnzzn — k_ _n zq 

K . U 
1 


and 


R ^m 


T (x)| = |_<2|cose| f / qJ (nkqsine)sln(^)cosa(q)dpdq (10-b) 

R. 0 
1 


n 


■3 3n‘ 


where is the peak chordwlse loading and a(q) is the geometric angle of 

attack. Equations (10) and (11) describe bounds on torque and lift or thrust 
noise, respectively, when they are used in eq. (2). 
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We have shown above that, if we are only allowed to deform the airfoil 
shape or the chordwise loading in such a way that the chordwise distance 
between points of equal thickness or equal loading is kept fixed, the 
maximum thickness and loading noise correspond to symmetrical positioning qf 
such points with respect to the same radius. Incidentally, in this case the 
thickness noise and loading noise are 90 degrees out of phase. 

To find the absolute maxima of thickness and loading noise, eqs. (9) and 
(10) will be used. It is assumed that the airfoil thickness and load distribu- 
tion functions are monotonic with respect to variable b. This assumption is 
satisfied in most cases of Interest. To be specific, thickness noise will be 
considered first. The right side of eq. (9) can be written as 


R m 

/ / nJ (nkn sine) sin dhdn 

T, n 2n 



= - / / 

R. 0 


nJ (nkn sine) sin dbdri 

n zri db 


( 11 ) 


where is the blade chord. To maximize the last integral, take 


j-=-A6(b-b) (12) 

db p 

where A is a function determined by the maximum thickness of the airfoil and 

b =b (n) is the value of b<b where sin(nb/2Ti) achieves its peak. If t(n) is 
p p o __ 

the thickness ratio of the airfoil, then A(n)=b^(n)t (n) . Using eqs. (9), (11) 
and (12) , we find that 

^ R ^ nb 

I'l' (x) I = 2 / nb t J (nknsine)sin(-y~2-) dq (13) 

n ' „ o n 2q 

max R . 


Similarly, to maximize the Integrals in eqs. (10-a) and (10-b), take 

= - B6(b-bp) (14) 

where B is a function determined by spanwise loading of the blade. If Cn (n) is 

2 2 

section lift coefficient of the blade, then B=p b n to C|,/2b . Equations 
(10-a) and (10-b) then give 
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|Q„S)I - , 

max R 


2 2 

p^cj R b^n nb 

/ — ^ J (nkrisine)sin(— ^)sina(n)dn 

Dp n n 


(14-a) 


l\(x)| = PqW I cose I / 

max 


R’ b ri c nb 

— ^ J^(nkr)s±ne)sin( 2 — ^)cosa (n) dn (14-b) 

R^ p " 


We have shovm that for all blades with a given thickness ratio t(n) at 
each radial station, the rms amplitude of the nth harmonic of thickness noise 
has the following bound 


^Trto|p^(x)|j< /2 p^n^(o^l4’^(x)| 

max 


(15) 


Similarly for all blades with a given section lift coefficient Cj^(n), the rms 
amplitude of the nth harmonic of- torque and thrust (lift) noise have the bounds 


47rro| p^(x) |j< /2 nk [Q^(x) 1 (16-a) 

max 

4irro|P]!^(x) |j<v^ nk|Tn(x)| (16-b) 

max 

respectively. 

Equations (12) and (14) show that the thickness function and chordwise 
load distribution function which generate maximum noise have infinite slope at 
the same two points which are symmetrically located with respect to the mid- 
chord. The corresponding thickness and chordwise load distributions are rec- 
tangular. Note also that jp^j is maximized by different thickness and load 

distributions for different n. In general, therefore the results of equations 
(15) and (16) are expected to be too pessimistic. 

Figure 2 shows some spanwlse aerodynamic data for a two-bladed helicopter 
rotor. These performance data were calculated by a strip theory-momentum 
analysis described in reference 4. The blade thickness ratio is 8 percent, 
the blade radius is 5 m and the chord is 0.4 m. The blade planform is rectangu- 
lar. The rotor rpm is 626. Figure 3 shows the calculated thickness and load- 
ing noise with the theoretical bounds obtained above. The chordwise load distri- 
bution at each radial position was obtained from the Garabedlan-Kom program 
(ref. 5). The airfoil section used is NACA 0008. The observer position is 50 m 
from the rotor center and 30° below the rotor (e=120°) . The bound for loading 
noise is obtained by adding the right sides of eqs. (16-a) and (16-b) . It is seen 
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that both bounds are very coarse although the bound on loading noise is not as 
pessimistic as that of thickness, noise. For the first harmonic level of loading 
noise, the reason for the theoretical bound being lower than the calculated level 
is not known- It may be due to the fact that the drag force (skin friction and 
wave drag) obtained from the Garabedlan-Korn program has a component normal to 
the chord which was used in the acoustic calcualtions. 

For B blades, substitute nB in all the equations derived above. 


THE EFFECT OF BLADE SWEEP 


We have shown that l'i'^(x)| is maximum when S^(n ,h)=constant. One way of 

reducing the level of the thickness and loading noise is blade sweep. This can 
be seen from eq. (8). The question arises whether a blade sweep can be select- 
ed which generates the least noise. We will show that among the blades with 
gradually Increasing sweep towards the tip, there is no optimal sweep. 

2 

To be specific, we take B^=-yri where y>0. The same argument holds as 

long as 38^/3ri < 0. From eq. (6), we have 

„ h 

R . 2/0 ™ 

’f (x) = / n J (nknsine) ^ dn / sin(-^) dh 

n R n Q 2n 

i 


= / g(n) dn (17) 

R. 

X 


where g(n)>0 is defined as 

h 

m 

g(n) = n J (nknsine) I 
n 


sin(“) dh 
-in 


(18) 


It is assumed, as before, that nb/2n£^ir. We note that as y Increases, so does 
the blade sweep. We have 


dl'l' (x)r R R , 

“ dy g(n)g(n' )cos[^(n -n' )]dndn’ 

R . R , 

1 1 


R n * 

= -ny / f (n^-n'^)g(n)g(n')sin[^(n^-n'^)]dndn' 
R. R. 

i X 


(19) 
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(20) 


If it Is assumed that sin ]>0, that is n^2Tr/y (R^-R^'^) , then 




For all practical angles of sweep, the above restriction on n is less strict 
than previously obtained restriction n_<2irn/b. The above result indicates that 
the levels of harmonics of the thickness noise decreases as the blade sweep in- 
creases. This result is also valid for loading noise. 

To test the validity of the above result, figure 4 shows the thickness 
noise spectra of three blades with increasing sweep. The tip Mach number is 
0.95 and the thickness ratios of all the blades are 8 percent. The observer 
is in the rotor plane and 50 m from the rotor center. It is seen that the above 
conclusion is indeed correct and should hold up to the 22nd harmonic. In fact it 
holds for much higher harmonics. The airfoil section used in the calculations 
is NACA 0008. 


CONCLUSIONS 


In this paper, bounds are established on thickness and loading noise of ro- 
tating blades. Only steady loading noise is considered which restricts the re- 
sults to high tip speeds. It is shown that only chordwise thickness and load 
distributions with midchord symmetry need be considered to establish these 
bounds. The resulting thickness and load distributions have infinite slopes at 
two points symmetrically located with respect to the midchord. Due to the fact 
that the amplitude of each harmonic of the spectrum is maximized, the resulting 
bounds are too coarse. A more appropriate approach may be to search for chord- 
wise thickness and load distributions which maximize overall acoustic level. 

It is also shown that sweeping the blade tips is beneficial in reducing 
the radiated noise. Also for blades with sweeps that increase towards the tip, 
there is no optimal sweep for minimum noise level. 
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(a) Thickness. 



Figure 1.- Chordwise distributions. 
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Figure 2.- Performance curves for two-bladed helicopter. 
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HARMONIC NUMBER 


(a) Thickness noise. 


(b) Loading noise. 


Figure 3.- Comparison of theoretical bound with calculated thickness 

and loading noise spectra. 
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